The impact of melt ponds on summertime microwave brightness temperatures and sea-ice concentrations

Sea-ice concentrations derived from satellite microwave brightness temperatures are less accurate during summer. In the Arctic Ocean the lack of accuracy is primarily caused by melt ponds, but also by changes in the properties of snow and the sea-ice surface itself. We investigate the sensitivity of eight sea-ice concentration retrieval algorithms to melt ponds by comparing sea-ice concentration with the melt-pond fraction. We derive gridded daily sea-ice concentrations from microwave brightness temperatures of summer 2009. We derive the daily fraction of melt ponds, open water between ice floes, and the ice-surface fraction from contemporary Moderate Resolution Spectroradiometer (MODIS) reflectance data. We only use grid cells where the MODIS sea ice concentration, which is the melt-pond fraction plus the ice-surface fraction, exceeds 90 %. For one group of algorithms, e.g., Bristol and Comiso bootstrap frequency mode (Bootstrap_f), sea-ice concentrations are linearly related to the MODIS melt-pond fraction quite clearly after June. For other algorithms, e.g., Near9OGHz and Comiso bootstrap polarization mode (Bootstrap_p), this relationship is weaker and develops later in summer. We attribute the variation of the sensitivity to the melt-pond fraction across the algorithms to a different sensitivity of the brightness temperatures to snow-property variations. We find an underestimation of the sea-ice concentration by between 14 % (Bootstrap_f) and 26 % (Bootstrap_p) for 100 % sea ice with a melt-pond fraction of 40 %. The underestimation reduces to 0 % for a melt pond fraction of 20 %. In presence of real open water between ice floes, the sea-ice concentration is overestimated by between 26 % (Bootstrap_f) and 14 % (Bootstrap_p) at 60 % sea-ice concentration and by 20 % across all algorithms at 80 % sea-ice concentration. None of the algorithms investigated performs best based on our investigation of data from summer 2009. We suggest that those algorithms which are more sensitive to melt ponds could be optimized more easily because the influence of unknown snow and sea-ice surface property variations is less pronounced

Optimal estimation of sea surface temperature from AMSR-E

The Optimal Estimation (OE) technique is developed within the European Space Agency Climate Change Initiative (ESA-CCI) to retrieve subskin Sea Surface Temperature (SST) from AQUA’s Advanced Microwave Scanning Radiometer—Earth Observing System (AMSR-E). A comprehensive matchup database with drifting buoy observations is used to develop and test the OE setup. It is shown that it is essential to update the first guess atmospheric and oceanic state variables and to perform several iterations to reach an optimal retrieval. The optimal number of iterations is typically three to four in the current setup. In addition, updating the forward model, using a multivariate regression model is shown to improve the capability of the forward model to reproduce the observations. The average sensitivity of the OE retrieval is 0.5 and shows a latitudinal dependency with smaller sensitivity for cold waters and larger sensitivity for warmer waters. The OE SSTs are evaluated against drifting buoy measurements during 2010. The results show an average difference of 0.02 K with a standard deviation of 0.47 K when considering the 64% matchups, where the simulated and observed brightness temperatures are most consistent. The corresponding mean uncertainty is estimated to 0.48 K including the in situ and sampling uncertainties. An independent validation against Argo observations from 2009 to 2011 shows an average difference of 0.01 K, a standard deviation of 0.50 K and a mean uncertainty of 0.47 K, when considering the best 62% of retrievals. The satellite versus in situ discrepancies are highest in the dynamic oceanic regions due to the large satellite footprint size and the associated sampling effects. Uncertainty estimates are available for all retrievals and have been validated to be accurate. They can thus be used to obtain very good retrieval results. In general, the results from the OE retrieval are very encouraging and demonstrate that passive microwave observations provide a valuable alternative to infrared satellite observations for retrieving SST

A new tracking algorithm for sea ice age distribution estimation

A new algorithm for estimating sea ice age (SIA) distribution based on the Eulerian advection scheme is presented. The advection scheme accounts for the observed divergence or convergence and freezing or melting of sea ice and predicts consequent generation or loss of new ice. The algorithm uses daily gridded sea ice drift and sea ice concentration products from the Ocean and Sea Ice Satellite Application Facility. The major advantage of the new algorithm is the ability to generate individual ice age fractions in each pixel of the output product or, in other words, to provide a frequency distribution of the ice age allowing to apply mean, median, weighted average or other statistical measures. Comparison with the National Snow and Ice Data Center SIA product revealed several improvements of the new SIA maps and time series. First, the application of the Eulerian scheme provides smooth distribution of the ice age parameters and prevents product undersampling which may occur when a Lagrangian tracking approach is used. Second, utilization of the new sea ice drift product void of artifacts from EUMETSAT OSI SAF resulted in more accurate and reliable spatial distribution of ice age fractions. Third, constraining SIA computations by the observed sea ice concentration expectedly led to considerable reduction of multi-year ice (MYI) fractions. MYI concentration is computed as a sum of all MYI fractions and compares well to the MYI products based on passive and active microwave and SAR products